Amplification of Curvature Perturbations in Cyclic Cosmology 
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We analytically and numerically show that through the cycles with nonsingular bounce the am- 
plitude of curvature perturbation on large scale will be amplified and the power spectrum will be 
redden. In some sense, this amplification will eventually destroy the homogeneity of background, 
which will lead to the ultimate end of cycles of global universe. We argue that for the model with 
increasing cycles, it might be possible that a fissiparous multiverse will emerge after one or several 
cycles, in which the cycles will continue only at corresponding local regions. 
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Recently, a cosmological cyclic scenario, in which the 
universe experiences the periodic sequence of contrac- 
tions and expansions |l'|, has been rewaked \2\, and 
brought the distinct insights into the origin of observ- 
able universe. There have been lots of studies on cyclic 
or oscillating universe models |i,|i,[l, 0,0,11,111, [13, 
[ll], 12], 13], also [14] for a review. In general, it is 
thought that the background of cyclic universe is homo- 
geneous cycle by cycle all along. However, it has been no- 
ticed that in the contracting phase the amplitude of cur- 
vature perturbation on super Hubble scale is increased, 
while is nearly constant in the expanding phase. Thus 
the net result is that the amplitude of perturbation is 
amplified, which seems to go along cycle by cycle [IQ] . 

Whether this amplification of perturbation actually oc- 
curs is interesting, since it will lead the global config- 
uration of cyclic universe unexpected colorful 161]. In 
this paper, we will analytically and numerically show 
the change of power spectrum of curvature perturbation 
through one or several cycles. We will simulate the effect 
of the increasing of its amplitude on the global configu- 
ration of cyclic universe, which might dramatically alter 
our conventional perspective for cyclic universe. 

We begin with the cyclic universe models in FiglH 
in which the nonsingular bounce is implemented by in- 
troducing a field with negative energy, which is mini- 
mally coupled to gravity. During the expansion and con- 
traction, the universe is dominated by a normal scalar 
field, which oscillates in a quadratic potential and offers 
matter-dominated background. With the shrinking of 
the scale factor, the ghost field becomes dominated and 
finally leads the bounce take place. In the inset of Fig.l, 
the corresponding state parameters uj = p/ p of models 
are plotted. It can be found that, due to the oscilla- 
tion of the dominated scalar field, w is oscillating with 
(w) ~ [l7|, which is matter-liked, during the contrac- 
tion and expansion. While near the bounce, uj crosses —1 
in a very short time. 

In general, the introduction of field with negative en- 
ergy is only the approximative simulation of a funda- 
mental theory below certain physical cutoff. Pointed by 
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FIG. 1: The evolution of scale factor in cyclic universe models. 
The blue line denotes that with increasing cycles. The red line 
denotes that with equal cycles. The time is scaled by ka which 
is an adjustable parameter to make the models be able to be 
consistent with the observable universe. 



[l8|, for a ghost with a minimal coupling to gravity, as 
same case in our paper, the cutoff A is constrained by 
observations of the diffuse gamma ray background with 
A < 3MeV. Therefore, we can estimate the validity 
of our model by consider the valid physical momentum 
space of perturbation. Given by the cutoff A, the up- 
per limit of the comoving wave numbers is fc^, ~ aK/kg 
where kc is in units of kg^ an adjustable parameter that 
can be determined by fit the time of expansion, which 
is 40/fcs, see Fig.l, to the age of observable universe, 
4Q/ks ^ I/Hq. When we use the present Hubble rate, 
Hq ^ lO^'^'^eV, we can find that, during the bounc- 
ing, the upper limit of the perturbations' comoving wave 
numbers is kc ~ 10^^ which is much larger than the 
largest wave number, 10^, which is chose in following 
figures', discussed in this paper. In this sense, the valid- 
ity of the background model used in this paper is well 
guaranteed. However, one should remember that the ap- 
pearance of phantom field is only an artificial approxima- 
tion. Recently, the nonsingular bounce has been obtained 
in nonlocal higher derivative theories of gravity [l9|,|20|. 
which can be ghostfree. 

The primordial perturbation generated during the con- 
traction with w ~ is scale invariant ^21], 22], [23]. There 
also are some studies how t he p erturbations go through 
such a nonsingular bounce [3, HI, [11], [13 ■ The cycle 



plotted as the blue line is increased, because there is an 
inflationary epoch after the bounce in each cycle, the cor- 
respondinK observable signals have been studied in e.g. 
[28], 29], 30]. The models with increasing cycle can be 
also obtained in e.g. 0j|3lj[liljOj '^^ which the entropy 
is increased. However, the results of change of power 
spectrum of curvature perturbation through cycles are 
not altered qualitatively by how the increasing cycle is 
implemented. The examples in Fig.l will only serve the 
purpose of numerical simulations in the following. 

We will regard the turnaround time as the beginning 
of a cycle, which can be denotes as i^ of the j*'' cycle. In 
each cycle the universe will orderly experience the con- 
traction, bounce, and expansion, and then arrive at the 
turnaround, which signals the end of a cycle. The begin- 
ning and the end of each phase of the j*^ cycle can be 



denoted as t^^. 



'•Ce' ''Ei' 



t^g, and the bounce is t^. 



The evolution of curvature perturbation under this 
cyclic background can be simply showed as follows. The 
motive equation of the curvature perturbation ^ in the 
momentum space is 



+ fc' ]uk = 0, 



(1) 



where Uk = zC.k |3l|,[32|, and [33| for details, the prime 
denotes the derivative for the conformal time 77 and z = 
-g.(^^)^/^. When fc^ :^ z" /z, i.e. the perturbations are 
deep inside the Hubble radius, it is obviously that Uk will 
osciUate with a constant amplitude. When fc^ ^ z"/z, 
i.e. the perturbations are on super Hubble scale, Eq.([T]) 
has general solution Uk, which gives Cfe — Ci + C2 J -p-, 
e.g. [33, where Ci and C2 are constant for fixed k. In 
general, for the contraction phase a ~ (ts—t)^, a ~ {jjb — 
T]) 1-" can be obtained. Thus during the contraction with 
n > 3 , as i -s> ts , the amplitude of C will be dominated 
by C2 term. 
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(2) 



since / I? 



J -^ ^ [rjB — 1]) 1"" , which is increased on 
large scale. 

While for the expanding phase a r^ t^, a ^ (±7]) 1-" 
can be obtained, where the plus is for g < n < 1 and 
the minus for n > 1. In these both cases the C2 term 
is decreasing. Thus the Ci term dominates Cfej which is 
constant for fixed k. This means that for a cycle of cyclic 
universe during the contraction C,k is increased on super 
horizon scale, up to the end of contracting phase in cor- 
responding cycle, while during the expansion it becomes 
constant. Thus the net result is that C,k on large scale is 
amplified, which is inevitable. 

The spectrum index of ^ obtained by Eq.([T]) is n^ — 1 = 



3- 



3n-l 



,e.g.[34 



inflation, and for 



which is scale invariant for n 3> 1, i.e. 
the contraction with w ~ 



I.e. 



During the contraction, when n > ^ the 
C2 term dominates, while when < n < :t, C is not 
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FIG. 2: The evolutions of perturbations through cycles given 
in Fig.l. The left upper and lower panels correspond to the 
models with increasing cycle and equal cycle in Fig.l, respec- 
tively. The black line denotes the evolutions of the Hubble 
radius ^j-mr, and the blue and red lines denote the evolutions 

of the perturbations with the modes 10~^fcs and fcs, respec- 
tively, which crosse the Hubble radius at different time. The 
evolutions of corresponding modes are given in the right up- 
per and lower panels, respectively. The evolutions of |$fcp 
are denoted with thick lines, while that of |("fcp are denoted 
with thin lines. We can see that C,k is increased cycle by cycle. 



dominated by the C2 term, i.e. its increasing mode, see 
Eq.(I2]). Thus in this case, it seems that the curvature 
perturbation is not amplified, n ~ 0+ corresponds to 
that in ^2], however, see [35| for the case with changed w. 
The net amplification for the perturbation modes, 
which are on super horizon scale in the f^ and j"* -I- 1 
cycles all along, can be estimated as follows. In the j*'* 
cycle, after the bounce, C, is given by Eq.®, which will 
be unchanged up to the end of the j cycle. Then the uni- 
verse enters into the j^^ + 1 cycle, during the contraction 
of the j*'' 4- 1 cycle, C, will continue to increase. Thus at 
certain time during the contraction, we have 



(^j+i(iJ+i) 



j+i _ p+i \ ^^^"^+1 



j+i _ ^i+i 
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'iJ + l„imV~'"^ + ^ 
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,B ''Ci , 
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l — 3nj 



exp 



(3n^+i - l)A/'^+i j^, 



1 - nJ+i 



Ce' 



(3) 



since H ^ ^/{tB — t), which means that for these modes 
still staying on super horizon scale the amplitude of the 
spectrum will be amplified with same rate after each cy- 
cle. In the third line 



AA^'+i = In 



^Ci ^Ci , 



(l-nJ+i)ln 






(4) 



is the e-folding number for the primordial perturbation 
generated during the contraction of the j*^ -I- 1 cycle. 
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FIG. 3: The left panel is the power spectrum of perturbations 
for model with increasing cycle in Fig.l at different time. The 
times chose to plot the power spectrum have been signaled in 
the figure, kst = 60 is that after the bounce in the ji"' cycle 
and kst = 160 is that after the bounce in the j"' -I- 1 cycle 
and k is in the units of ks as in Fig.l. The right panel is 
that of perturbations for model with equal cycle in Fig.l. We 
can see that after one cycle the spectrum will be redden on 
corresponding scale. 



When w ~ 0, i.e. n^+i ~ |, Eq.© becomes C^+^{P+'^) ~ 

e^ TJpg, which is consistent with that in |15l.llq. 

We numerically show the evolutions of the perturba- 
tion modes Cfc in details under the backgrounds in FigHJ 
which are plotted in proper time. This actually can be 
manipulated by numerically solving the equation of the 
metric perturbation $fc with cj. = 1. Then the evolu- 
tion of Cfc is obtained by using C = f( ^'i+t^'^ ) + *• We 
can clearly see that (^k is increasing during the contrac- 
tion and is constant during the expansion, which thus 
is amplified cycle by cycle. The results are consistent 
with above discussions, and also the numerical results in 
[ill.d^l for one cycle. 

We neglected the effect of entropy perturbation in this 
simulation. However, the inclusion of the entropy per- 
turbation will not change the result essentially, since it is 
generally not important. Though in the background evo- 
lution of Fig.l, the nonsingular bounce is implemented by 
introducing the field with negative energy, the numerical 
results obtained might be applicable for other cyclic mod- 
els with nonsingular bounce. The differences are that in 
these models around the bounce c^ might change with 
time and the higher order terms of $ appear. However, 
if the correction for the equation of $ is important only 
around the bounce, the results are not expected to be 
altered by these difference. 

The power spectrum of C is 7^c(^) — ^iCfcP- We 
assume that the power spectrum of ( after the bounce 
of the j*^ cycle is 'P(;{k,t^g). Thus in term of Eq.©, at 
certain time P~^^ during the contraction of the j*'' -1- 1 
cycle, the power spectrum of the perturbation modes, 
which are all along on super horizon scale, is given by 
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which means that the amplitudes of these modes are am- 
plified, but the spectrum shape is unchanged. There are 



also the perturbation modes, which will enter into the 
horizon during the j*^ cycle and then leave it during the 
contraction of the j*'* + 1 cycle, see the red lines in right 
panels of Fig. 2. Inside the Hubble radius, Uk is constant. 
Thus we have Cfc — Uk/z '^ Uk/a. This means that the 
amplitude of Cfe is decreasing during the expansion and 
increasing during the contraction on subhorizon scale. 
The power spectrum of these modes can be estimated as 



V'+\k,P+') 



^^^ . ~ ^g' > t^ + ' - P- 

(P _ fJ ^2r^. V ^ 



(t'Ee-t's) 






kg k J 
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(6) 



W'-k) 



for t-jjl < P^^ < t'J'^ , where A;^, kg are the modes 
crossing the horizon at the bounce time t'g. In general, 
the change of spectrum is not exactly power law. The 
second line is obtained only when k ^ kg, ^g . We 
can see that for these perturbations modes, not only the 
amplitude of the spectrum is increasing after each cycle, 
but also the shape of the spectrum is changed with 



An^ 



An'+^ - 2 



2nij 



,i+i 



- 1 



(7) 



where 



and 



.J+i 



denote that in the expanding phase 



of the i*^ cycle and that in the contracting phase of 
the i*^ + 1 cycle, respectively, which are generally not 
equal. In general, for the interesting value of n^ and 



3 < •ni'fj < 1, AriQ is negative, thus the spectrum will 
redshift. However, it is possible that An^ ~ by choos- 
ing special n^ and n^ , which might has interesting 
application, e.g.[3g. In this calculation, we neglected 
the effect of the transfer function, which is dependent on 
the matter content in corresponding cycle. However, the 
evolutive behaviors of spectrum obtained here are not al- 
tered qualitatively. When n'^ = n;^ = n, the shape of 
spectrum is changed as An^ ~ ^^Ii ; which is consistent 
with the result in j37| . 

We numerically show the change of the power spec- 
trum through one cycle in Figl3] under the backgrounds 
in FigdJ The power spectrum VUkjtg) of ( after the 

bounce in the j cycle is scale invariant, since n ~ ^ . Dur- 
ing the contraction of the j*^ + l cycle, for the model with 
equal cycle the spectrum will redshift. While for that 
with increased cycle, the shape of spectrum on large scale 
is unchanged, i.e. still scale invariant, only its amplitude 
is amplified, since these modes are on super horizon scale 
during the contraction of the j*^ cycle and j*'' -I- 1 cycle 
all along, however, the spectrum of the modes on middle 
scale will redshift, and the rate of tilt can be estimated as 
An^ ~ 6. The spectrum on small scale is scale invariant, 
because these modes are newly generated in the j*'* -I- 1 
cycle. These results are consistent with Eqs.® and ^, 
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FIG. 4: C,{x) in position space, which reflects the inhomo- 
geneity of background. The left panel is for the model with 
increasing cycles at the time kst = 260, while the right panel 
is for that with equal cycles at the time kst — 110, and the 
length scale is in unit of the Hubble radius at that time. 



which obviously reflect the changes of power spectrum 
through cycles. 



1 /n 

Therefore, for -i < n < 1 the amplitude V^ of per- 



1/2 
C 



1 



turbation on large scale will inevitably arrive at V 

at certain time of cycles. We only show it '^ 0.1 in Fig|3l 

since the linear perturbation approximation is not reli- 



able when V, 



1/2 

c 



1, and in this case, the nonlinear effect 



will become important and the coupling between modes 
has to be considered. However, it can be expected that 

1/2 

the enhancement of nonlinear effect will bring V/ to 1 
faster. This means that on the corresponding scale the 
perturbation will be large enough so that will destroy the 
homogeneity of background. We illustrate this effect of 
the perturbation on background by transforming V/ in 
Fig. 3 into ({x) in position space plotted in Fig. 4. 

1/2 

We can see that when V/ -^ I the universe is frag- 
mentized into large number of small regions. In this case, 
the initially global homogeneous universe will become 
highly inhomogeneous. Thus it can be hardly imagined 
that the different regions of global universe will evolve 
synchronously, even if it is homogenous in the previous 
cycle. Thus the cycle of global universe will ultim ately 
end. The similar phenomenon was also discussed in [4ll |. 
however, in which since the bounce is singular the special 
mode mixing has to be applied for the amplification of 
C: e.g.[33|,[30|,[4O]. Here, the bounce builded is nonsin- 
gular, thus the evolution of perturbations through cycles 
can be exactly studied by numerical method. 

In general, as long as the curvature perturbation on 



large scale is amplified cycle by cycle actually occurs, it is 
hardly possible that the global universe is eternal cyclic. 
We only consider the evolution of curvature perturba- 
tion, the effect of the increasing of other perturbations 
on cyclic universe has been also studied, e.g.[43|- These 
results mean that the reliability of some models of cyclic 
universe might need to be reevaluated. 

However, it can be noticed that for the universe with 
increasing cycle, the regions split can be larger. When 
the length of local regions split is larger than the Hub- 
ble scale at the corresponding time, these regions will 
possibly evolve independently, as long as inside the cor- 
responding regions the background is homogeneous. In 
principle, such different local regions correspond to dif- 
ferent universes [4j] , [4Jj , each of which is controlled by 
local physical equations and might be fragmentized again 
after itself cycles. In this case, the number of local uni- 
verses will increase cycle by cycle. Thus we can have a 
cyclic multiverse scenario, as argued in [i5|,[la]. This ar- 
gument in some sense indicates that in cyclic cosmology 
a gradually increasing cycle is significant for the contin- 
uance of cycle, however, in this case, the cycle will con- 
tinue only at local regions, which is homogeneous. This 
scenario can be distinguished from that in chaotic eter- 
nal inflation [45 . 46] , in which the inflationary multiverse 
is induced by the large quantum fluctuation of inflaton 
field, which occurs efold by efold. However, the cyclic 
multiverse is induced by the cyclic amplification of per- 
turbation on large scale, which is in classical sense, and 
occurs cycle by cycle. 

In conclusion, we have analytically and numerically 
showed that through the cycles with nonsingular bounce 
the amplitude of curvature perturbations on large scale 
will be amplified and the power spectrum will be redden. 
In some sense, this amplification will eventually destroy 
the homogeneity of background, which might lead to the 
ultimate end of cycles of global universe. However, it 
can be argued that for the model with increasing cycle, 
the global universe will possibly evolve into a fissiparous 
multiverses after one or several cycles, in which the cycles 
will continue only at corresponding local regions, inside 
which the background is homogeneous. 
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